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Abstract 

Here we describe the quantum limit to estimation of the spacetime 
metric, or equivalcntly the quantum hmit to measuring the classical 
gravitational field. Specifically, wc write down the optimal quantum 
Cramer-Rao lower bound, for any single parameter describing a met- 
ric for spacetime. Four key examples are demonstrated covering a 
broad range of relativistic phenomena. Wc describe quantum limited 
estimation of the mass of a black-hole, the acceleration of a uniformly 
accelerating observer, the amplitude of a gravitational wave and the 
expansion parameter in a cosmological model. The standard time- 
energy uncertainty relation and the Heisenbcrg uncertainty relation 
are special cases of the uncertainty relation for the spacetime metric. 
The uncertainty relation takes a particularly simple and revealing form 
when the measurement region is made sufficiently small. We use the 
locally covariant formulation of quantum field theory in curved space- 
time, which allows for a manifestly spacetime independent derivation. 
The result is an uncertainty relation applicable to all causal spacetime 
manifolds. 

1 Introduction 

The geometry of spacetime is determined by physical measurements made 
with clocks and rulers or, more generally, using quantum fields, sources and 
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detectors. As these are physical systems, the ultimate accuracy achievable 
is determined by quantum mechanics. In this paper we obtain a parameter- 
based quantum uncertainty relation, bounding knowledge of the spacetime 
metric. 

The gravitational field is treated entirely in accordance with classical gen- 
eral relativity. In general relativity, the gravitational field is a manifestation 
of the geometry of spacetime, which in turn is described by a metric. The 
metric essentially gives the proper times and proper distances between events 
in spacetime. The infinitesimal proper distance between two spacetime events 
is given by pQ: 

ds"^ = g^y{x)dx'^dx^ (1) 

Where g^u{x) is the metric tensor, with indices /x, z/ = 0, 1, 2, 3 for the time 
and three spatial components. The dx^ are the infinitesimal coordinate dis- 
tances. We assume the Einstein summation convention, where repeated up- 
per and lower indices, are to be summed over, for example in equation (1). 

We know from Heisenberg that there is an uncertainty relation between 
position and momentum. It states that the product of the uncertainty in 
position and the uncertainty in momentum, must always be greater than a 
constant. We write down Heisenberg's uncertainty relation, in terms of the 
variances of position and momentum: 

((Ax)^)((Ap)2) > ^ 

This relationship bounds knowledge of the position and momentum of a 
quantum system, where h is given by the reduced Plank's constant. The 
Heisenberg uncertainty relation is derived in standard quantum mechanics, 
where position and momentum are both represented as Hermitian operators. 
A similar relation also exists between time and energy: 

((5tf)((A^f)>| 

Unlike position, time in standard quantum mechanics is given by a classical 
parameter t with classical variance The time-energy uncertainty 

relation can be considered as an example of quantum parameter estimation. 
There, one tries to estimate a classical parameter, in this case time, by making 
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measurements on a quantum system, which depends on the parameter of 
interest. 

The most common way to make quantum mechanics compatible with 
classical relativity, is to demote position to a parameter, just like time is in 
the previous example. Physical systems can then be thought of as living on, 
and interacting with, the classical spacetime manifold. In this relativistic 
context, the Heisenberg uncertainty relation should also be framed in the 
parameter estimation context. 

This way of thinking was used by Braunstein, Caves and Milburn [2], 
to develop optimal quantum estimation for spacetime displacements, in flat 
Minkowski spacetime. In spacetime, not only can one move a fixed proper 
distance or time, one can also boost and rotate. Quantum parameter esti- 
mation was thus also developed for the parameters corresponding to these 
actions |2]. The results were developed with the quantised electromagnetic 
field as the measurement system. They show that the estimate of spacetime 
translations may be made more accurate, if the uncertainty in the number 
operator is made very large. 

Before we describe quantum parameter estimation in more detail, we ask 
the question: can any insight be found by simply applying the Heisenberg 
uncertainty relation, in its parameter based form, directly to a proper dis- 
tance? It was in this way that Unruh [3] derived an uncertainty principle 
for the Qii component of the metric tensor. The calculation was simple and 
made the important point that in general relativity the coordinate system 
is arbitrary and so there should only be quantum uncertainty in the proper 
time and the proper distance. As these are in turn related to the metric via 
equation (1), any uncertainty in the proper distance is equivalent to uncer- 
tainty in the metric. By applying the Heisenberg uncertainty relation to a 
proper distance, a simple but insightful uncertainty relation was found for 
one component of the metric. We write down the Unruh uncertainty relation 
for the XX component of the metric, in terms of variances as: 

((5^n)^)((Ar")^) > ^ 

The key finding of this uncertainty relation is the inverse proportionality 
to V"^, the square of the four-volume of the measurement. The conjugate 
variable to gu is the corresponding component of the quantised stress-energy 
tensor T^^, in this case the pressure in the x-direction. 
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In this paper we consider a more general context for deriving an uncer- 
tainty relation for the metric, by formulating it as a problem in quantum 
estimation theory. The metric Qfiuix) is defined for each point x on the man- 
ifold. If our quantum measurement system occupies some four-volume V, 
then the system will depend on the the metic at every point in that region. 
If we consider the metric to be an arbitrary function on the manifold, then 
we need to estimate an infinite number of parameters to define it completely. 
Instead we will consider regions of spacetime, which can be described by 
metrics, defined by a finite number of parameters 61,62 ■■■ On- For exam- 
ple the Schwarzschild metric, which describes the spacetime around a static 
non-rotating blackhole, is defined by only one parameter M, the mass of the 
blackhole. The task is to estimate the parameters by making measurements 
on physical systems living on the spacetime manifold. 

The general schema for quantum parameter estimation is represented in 
figure [1] An initial quantum state, represented by a density operator pq, is 
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Figure 1: The scheme for quantum parameter estimation. 

evolved by a unitary transformation U (6) dependant on the parameter 6 of 
interest. Measurements are then made on the system, producing measure- 
ment results u, which are in turn fed into an estimator 6{u) of the parameter. 

We consider generalised measurements, known as positive operator valued 
measures. These are Hermitian operators with the property: 

j E{uj)duj = 1 

Where u is the measurement result and 1 is the identity operator. The 
outcomes of a particular measurement will follow a probability distribution 
p{uj\Q) conditional on the actual value of the parameter G. The probability 
distribution for the outcomes u can be calculated as: 

p(u;|0) = Tr(p(0)EH) (2) 
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Where p(0) = U{Q)poU^ (Q) is the state after the 9 dependant interaction. 

The problem of estimating the parameter 6 is essentially that of choosing 
a value 9 to fit the probability distribution p{uj\9) to the actual measured 
values. A common example is the so-called, maximum likelihood estimator, 
which is the choice of 9 which retrospectively maximises the probability, of 
the observed measurement values. 

The variance of any estimate of the parameter 9, based on the distribu- 
tion of the observed measurements, will be bounded below by the classical 
Cramer- Rao lower-bound [2]: 



F(0) 

where F{Q) is the Fisher information for the measurement, given by: 

'dp{uj\9) 



F{Q) = / du 



piu\Q) \ 39 



=e 



The key ingredient is the rate of change of the probability distribution with 
respect to the parameter. A larger response due to a change in the parameter 
naturally gives rise to better estimation of the parameter. The rate of change 
of the probability distribution can be calculated as: 



89 

Where h is the generator of the unitary transformation U{9). We have de- 
scribed quantum parameter estimation in the Schrodinger picture where the 
density operator depends on the parameter. Due to cyclicality of the trace 
used to calculate both the probability distribution and its rate of change, 
the Heisenberg picture can equally be used, which can be quite useful in the 
relativistic setting. The key component is the generator h which generates 
changes in the probability distributions due to changes in the parameter. 
This can be seen directly for pure states, when one optimises over all pos- 
sible measurements. The optimal quantum Fisher information is then given 
by: 
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Where ((A/i[0])^) is the quantum variance of the generator h taken in the 
initial state. The optimal quantum Cramer-Rao lower bound can then be 
expressed as: 

One can now easily construct examples by simply identifying parameters 
and their corresponding generators. For example, the Hamiltonian is the 
generator of time translations, which gives the uncertainty relation presented 
in the introduction. Another important example is the number operator and 
phase which is the basis of Heisenberg limited phase estimation. 

2 An Uncertainty Relation for the metric g^^, 

The physical systems we consider here are quantum fields, for example the 
Dirac field for electrons. For clarity we restrict to the free scalar field. The 
generalisation to interacting and higher spin fields, should follow a similar 
argument to the one we present here for the scalar field. We consider only 
measurements on scales large enough, such that the quantisation of the grav- 
itational field itself can be safely ignored. We also work in the test-field 
approximation, where the gravitational field of the probe is ignored com- 
pletely. For example one might choose to ignore the gravitational field of 
the laser in a gravitational wave interferometer. It is important to note how- 
ever, that the action of the probe fields on the spacetime manifold will only 
strengthen the bound we derive below. That is, if a high energy field were 
used to probe the structure of spacetime, then the back-reaction of the field 
on the spacetime, would further limit the accuracy of the measurements, pre- 
venting our inequalities from being saturated. The problem of back-reaction 
when measuring the structure of spacetime has already been studied in some 
detail 0]. 

Apart from the test-field approximation, the quantum theory for the fields 
we consider here takes classical general relativity fully into account. We em- 
ploy the locally covariant formulation of quantum field theory on curved 
spacetime [5]. Local covariance refers to the global property, of having the 
same physics under all coordinate systems, i.e. general covariance, together 
with the local property that if two localised spacetime regions are equivalent, 
then they should describe the same physics regardless of whats happening at 
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distant regions of the universe. In this sense the locally covariant formula- 
tion fully takes into account the locality and covariance of general relativity. 
We only omit causality-violating spacetimes, which removes certain causal 
pathologies associated with closed time-like curves and time machines. Our 
results will hold for all causal spacetime manifolds. 

A key result of the locally covariant approach, is the calculation of how 
quantum observables respond to changes in the metric. Say we believe some 
particular region of the universe to be well described by a metric g^viO) 
depending on N parameters 6 = (^1,^2, • • • ,0n)- For this case the locally 
covariant approach can be used to calculate how any observable E{6) will 
respond to a change in any one of the parameters. The response is evaluated 
as the rate of change of the observable with respect to the parameter and 
as we noted in the introduction, this is just what we need to calculate the 
quantum Cramer-Rao lower bound. 

We first consider an arbitrary region of spacetime and divide it into three 
sub-regions as shown in figure 1. The lower region labeled is the region 
where the field is prepared. It is in the causal past of the upper region 
labeled where measurements are performed. The intermediate region 
Ml is inside the intersection of the causal past of A^"*" and the causal future 

We consider a diffeomorphism 0^ which smoothly deforms the metric 
inside the region A4 and acts like the identity everywhere else. The effect 
of this change on any observable E{6) in the measurement region can be 
calculated by the action of a map, known as the relative Cauchy evolution 
[5]. First we construct the observable E{6) for a particular value of the 
parameter, 6 = Q, and then use the relative Cauchy evolution to calculate 
the rate of change with respect to 6. Specifically the diffeomorphism chosen 
acts on the metric in the region of interest as: 

The number G can be thought of as the actual or expected value of the 
parameter. If we have an unbiased estimator 6 of the parameter 6 then the 
expected value of the estimate (9) = Q. 

We construct the measurement operators out of polynomials of localised 
field operators. For a particular measurement E{u\9), conditional on a pa- 
rameter 6 in the metric, we can construct the probability distribution of the 
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Figure 2: Spacetime diagram. The middle region is described by a metric 
QiiuiO) and it is here that one wishes to obtain the value of 6. 



p{uj\d) = Tt(pE{cu\9) 



outcomes u as: 



where p is the density operator representing the state of the system. The rel- 
ative Cauchy evolution gives the rate of change of the probability distribution 
with respect to the parameter as 0: 



dp{uj\e) 



d9 



TTUp[E{co\Q),Pm 



(3) 



where the operator P{6) is given by: 



P{9) 



dpT- 



M 



de 
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Here T^'^ is the renormalised stress-energy tensor of the scalar field used by 
the measuring device. Ai is the region of interest and dfi is the volume-form 
induced by gn^{Q). We identify the operator P{6) as the infinitesimal genera- 
tor of changes in the probability distribution due to changes in the parameter 
9. The cyclicality of the trace and the simple commutator form of the ex- 
pression in relation (1) means that the dependance can either be thought 
of as in the observable or the state. The standard techniques of quantum 
parameter estimation, as described in the introduction, can be immediately 
applied. For pure states, optimisation over all possible measurements, yields 
the following uncertainty relation: 

{iser){{AP[Q]r) > I (4) 

This relation gives the optimal quantum limit to measuring any parameter 
of a spacetime metric. It is locally and covariantly defined and applies to all 
causal spacetime manifolds. To demonstrate the ease with which it can be 
applied, we shall now present a number of examples. 



3 Estimating the Mass of a Black Hole 

As our first example, we shall consider the problem of measuring the grav- 
itational field, outside a spherically symmetric, non-rotating, massive body. 
The metric for the empty space outside such an object, is given by the 
Schwarzschild solution to Einstein's field equations [I]. In Schwarzschild 
coordinates, with units such that the speed of light c = 1, the metric be- 
comes: 



ds^ = -[1 — \dt^ + ll — \ dr^ + r^de^ + sin^ Odcj? (5) 



where M is the mass as observed from infinity and G is the gravitational 
constant. These coordinates have the intuitive features that; surfaces of 
constant t and r have area given by 47rr^ and the metric (5) becomes the 
metric of an inertial observer in fiat spacetime, as r becomes large. We choose 
the region of interest M. to have duration Lt = t2 — ti, length Lr = r2 — ri 
and solid angle Lq = J^^ J^^ sm6d6d(f) in Schwarzschild coordinates. 

The uncertainty relation (4) is defined in terms of a four-volume integral 
of a scalar quantity. Due to this, it is invariant under general coordinate 
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transformations. We shall take full advantage of this invariance and change 
to more well behaved coordinates, in order to evaluate the uncertainty re- 
lation. We transform to outgoing Eddington-Finkelstein coordinates, which 
are appropriate for describing photons propagating out of the gravitational 
potential pQ. The outgoing null geodesies are labeled with a new coordinate 
U = t — r* where: 

r* = r + 2MG In \r/2MG - 1 1 

In the coordinates {U, r, 6, (p) the metric takes the outgoing Eddington-Finkelstein 
form [Ij: 

(2MG\ 
1 j dU^ - 2dUdr + r\de^ + sin^ Odct)'') 

For this example we shall choose a state, such that the variance of the stress- 
energy tensor, is constant throughout the measurement region, when ex- 
pressed in Eddington-Finkelstein coordinates. In a different state the vari- 
ance might not be constant in these coordinates. The behaviour of the vari- 
ance, over the measurement region, will affect the explicit dependence of the 
uncertainty relation, on the shape and size of the measurement region. With 
the choice of state and measurement region we have made, the uncertainty 
relation (4) becomes: 



((5M)^)((Ar"")^) > 



h/c 



GLtLn[rl-rl 



where we have reinserted the correct factors of c. We see that the Fisher 
information depends on the variance of the energy-density, in Eddington- 
Finkelstein coordinates. The Fisher information also increases, as one in- 
creases the size of the measurement region. This makes for an increased 
effect of the mass M and hence a better estimate. The dependance on the 
distance from the mass is clearly stronger, than that of time or solid angle. 
This is due to the fact that the gravitation field varies in this direction. One 
achieves better estimation therefore, by extending the apparatus out in this 
direction. By defining our measurement region in Schwarzschild coordinates, 
we have chosen a particular shape for the apparatus. Although the Fisher in- 
formation does not depend on the coordinates, it will depend on the shape of 
the region. As we shall see, the underlying dependance on the four- volume of 
the measurement, will be the inverse dependence seen in the Unruh relation 
(3). 
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4 Estimation of Proper Acceleration 



We now give an example, for the use of the uncertainty relation (4) in flat 
spacetime. In this example it is the property of an observer we are inter- 
ested in, not the spacetime itself and so the observer effectively fixes the 
coordinate system. This example is perhaps best interpreted as a constant 
acceleration drive spacecraft, trying to determine it's own acceleration by 
local measurements. The natural coordinate system, for a uniformly acceler- 
ating observer, is know as the Fermi- Walker transported orthonormal tetrad. 
It is the coordinate system carried by the accelerating observer, in which he 
is instantaneously at rest. This coordinate system has limits, it can only be 
extended at most a distance from the observer, where a is the proper 
acceleration. The fiat Minkowski metric in these coordinates can be written 
as tU: 

ds' = -{1 + aende? + {der + {de? + {der 

Here we take the region of interest to be a 4-cube, centred at the origin, in 
these coordinates with sides of length L^i, L^2, L^s and duration L^o. Once 
again we choose a state where the variance of the stress-energy tensor is con- 
stant, in the observers coordinates. The uncertainty relation for the proper 
acceleration a is then given by: 

We see that better estimation is achieved by making the apparatus larger and 
having a bigger acceleration. However, if we make the apparatus as large as 
possible, by inserting the constraint < ^ we find: 

{{S~amAf'y)>(-^^y (7) 

Now the inverse is true, if one has reached an acceleration, such that the 
size of the apparatus has been limited, then larger accelerations will produce 
worse estimation. This can be seen in the first instance, from the stronger 
dependence on size as compared to acceleration. 
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5 Quantum Limited Gravitational Wave De- 
tection 



We now consider another example of the uncertainty relation (4), in estimat- 
ing the amplitude of a gravitational wave. In the linear approximation of 
Einstein's equations, the metric can be written as the flat Minkowski metric 
plus a small perturbation. 
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The simplest solutions to the linearised Einstein equations are plane-waves. 
We write down these wave solutions, in the transverse traceless gauge, as [1] : 

V = ^[A^.expikax")] (8) 

The amplitude of the wave A^iy has only two independent components cor- 
responding to the two polarisation states. For our example here, we shall 
consider a gravitational wave propagating in the z direction and linearly po- 
larised along the diagonal directions in the x-y plane. For this case the only 
non-zero components of the metric (8) are [1]: 

As the wave is a small perturbation Ay. ^ 1 we approximate the four- volume 
element, as being the volume element of the flat metric. For our measurement 
region we take a 4-cube centred at the origin, in Minkowski coordinates, with 
length LxiLy, Lz and duration Lj. Again we choose a state where the variance 
of the stress-energy tensor is constant. When estimating the amplitude the 
uncertainty relation then becomes: 

{{5A^f){{^f^yf) > ^'"^^ 



GALlLlsinY^ooLt] sin^f^wL,] 

One can see that larger frequency waves will produce worse estimates for 
the amplitude. Also the Fisher information oscillates to zero whenever the 
transverse length or duration of the apparatus is a multiple of the wavelength. 
Consider a ring of atoms sitting in the plane perpendicular to the direction 
of propagation of the wave. Once a complete wavelength has passed they will 
distort and then return to their original position. If a measurement is only 
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made at this point then no effect will be observed, and hence no information 
on the amplitude will have been received. A limitation to the uncertainty 
relation (4) in dynamical situations like this one, is that it is based on a single 
measurement. Generalising it to take into account continuous measurements 
is then the logical next step. It will require analysing the Fisher information 
matrix and will result in a spectral uncertainty relation [B]. 



6 Estimating the expansion of the universe 

For our final example we shall consider the spatially flat Friedmann- Robertson- 
Walker cosmology. This is a universe filled with a uniform density of galaxies. 
At any instant in time, in the co-moving frame of the galaxies, the universe 
looks the same everywhere (homogeneous) and in all directions (isotropic). 
The metric for this universe is given by 

ds^ = -df + a^{t) [dx'^ + dy"^ + dz^] 

where t is the proper time of an observer co-moving with any of the galaxies. 
The spatial coordinates x, z describe the homogeneous and isotropic sur- 
faces of constant proper time t. The function a(t), known as the expansion 
parameter, is the ratio of the proper distance between any two galaxies at 
the initial time t = and the time t. 

During an infinitesimal duration of proper time dt a photon will travel 
the distance drj = It is convenient to use this as the time parameter. 
We shall consider a universe dominated by matter, in which case the metric 
becomes; 

ds" = "^{l- cos riY[-dri^ + dx^ + dy'' + dz^] 

where rj runs between at the beginning of expansion to 27r at the end of 
recontraction. We wish to estimate the parameter flmax which controls the 
maximum size the universe reaches before recontraction commences. The 
uncertainty relation for this parameter becomes: 
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The integral can be performed analytically but perhaps becomes somewhat 
clearer if expressed in terms of the density of mass-energy p: 

^{{^amax) )((Af ) ) _ ^3^^^2.j^2>^^jna^ JJ?2^^^_2^^^) j 

Better estimation is made by running the experiment for longer time 77. How- 
ever for fixed duration, better estimation is made during the period of low 
energy-density, which occurs at the point of maximum expansion. 

7 Estimation of Proper Time and Proper Dis- 
tance 

Consider now making measurements in a local inertial frame, sufficiently 
small such that the metric can be approximated as constant, over the mea- 
surement region. The only parameters to estimate are then the components 
of the metric. For the first component of the metric, the uncertainty relation 
(9) becomes: 

{{SgoomA [ d,if''f)>re (9) 

JM 

The integral of the energy density over the spatial component of the four 
volume is the Hamiltonian: 

/ ci//rO° = ,/{g^) I dtH{t) 

and so the uncertainty relation (9) becomes: 

((5^00)^) ((A j dtH{t)f) > ^ (10) 

Now consider the proper time of an observer stationary in this inertial coor- 
dinate system. The proper time in this case is given by: 
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The coordinates have been chosen so there is no uncertainty in the coordinate 
time. Using the delta method, the uncertainty in the metric is then related 
to the uncertainty in the proper time by: 

moor)^M9oo){{srr)/t' 

Inserting this relation into equation (10) gives the time energy uncertainty 
relation generalised to curved spacetime in a local Lorentz frame: 

i(6rmAi j dtmf) > (11) 

If we assume a time independent Hamiltonian and specialise to flat spacetime 
the uncertainty relation (11) reduces to: 

{{STmAHf) > ^ (12) 

This demonstrates that the standard time energy uncertainty relation is a 
special case of the metric uncertainty relation (4) . Using a similar argument 
on can also derive an uncertainty relation for proper distance X. 

{{sxr){{Apr) > ^ (13) 

This is the parametric version of the Heisenberg uncertainty relation where 
P is the momentum in the direction of the displacement. 



8 Simplified Metric Uncertainty Relation 

Let us consider the uncertainty relation (9) for an arbitrary component of 
the metric: 

((%.)^)((A / dfifn')>n' (14) 

Jm 

Where in this case there is no sum over the repeated indices. We further 
simplify the uncertainty relation (14) by assuming the stress-energy tensor 
is also constant over the measurement region. In this case one is left with a 
particularly simple form for the metric uncertainty relation given by: 
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where V is the four- volume of the measurement. This form of the fundamen- 
tal dependance confirms the earlier relation (3) found by Unruh. Indeed the 
same restrictions were needed to derive (3) as were used to produce (15). 

9 Conclusion 

In this paper we have presented the optimal quantum Cramer-Rao lower 
bound for parameters describing a metric for spacetime. Our specific deriva- 
tion applies for pure states of the scalar field on an arbitrary causal spacetime 
manifold. We give four important examples covering the full gamut of rel- 
ativistic phenomena. We described quantum estimation for; the mass of a 
black-hole, the acceleration of a uniformly accelerating observer, the ampli- 
tude of a gravitational wave and the expansion parameter in a cosmological 
model. In all these examples the fundamental dependance can be seen as 
the inverse proportionality to the four- volume of the measurement. This 
dependance can be seen explicitly when one makes the measurement region 
sufficiently small, in agreement with earlier work on the subject |3]. 

The methods developed here can easily be applied to many situations in- 
volving the measurement of gravity. The uncertainty relation (4) can be used 
to benchmark the optimality of experimental proposals involving high pre- 
cision measurements of gravitational phenomena. By evaluating the form of 
the uncertainty relation (4) for particular scenarios, one can ascertain which 
parameters have the biggest impact on the measurement results. Hence one 
can better understand the tradeoff between cost and impact when attempting 
to optimise experimental design. 
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